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Abstract 

The analytical vectorial structure of non-paraxial four-petal Gaussian beams(FPGBs) in the far 
field has been studied based on vector angular spectrum method and stationary phase method. 
In terms of analytical electromagnetic representations of the TE and TM terms, the energy flux 
distributions of the TE term, the TM term, and the whole beam are derived in the far field, 
respectively. According to our investigation, the FPGBs can evolve into a number of small petals 
in the far field. The number of the petals is determined by the order of input beam. The physical 
pictures of the FPGBs are well illustrated from the vectorial structure, which is beneficial to 
strengthen the understanding of vectorial properties of the FPGBs. 
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I. INTRODUCTION 



Recently, a new form of laser beams called four-petal Gaussian beams(FPGBs) has been 
introduced and its properties of passing through a paraxial ABCD optical system have been 
studied Subsequently, the propagation properties of the FPGBs have attracted consider- 
able interest due to its potential applications. Gao and Lii studied its vectorial non-paraxial 
propagation in free space based on vectorial Rayleigh-Sommerfeld diffraction integral for- 
mulas in 2006 In 2007, the propagation of four-petal Gaussian beams in turbulent 
atmosphere was investigated by Chu et al jsl- In 2008, Tang et al explored diffraction 
properties of four-petal Gaussian beams in uniaxially anisotropic crystal by virtue of the 
paraxially vectorial theory of beam propagation {4]. Yang et al reported the propa^tion 
of four-petal Gaussian beams in strongly nonlocal nonlinear media in the next year [^]. In 
the meanwhile, the vectorial structure of lots of beam with different patterns and polarized 
status is illustrated in the far field by means of vector angular spectrum method J6|, which is 
a useful tool to resolve the Maxwell's equations, and stationary phase method [7|, Isj, which 
uses the asymptotic approximation approaching some kind of difficult integral. Based on 
these two methods mentioned above, Zhou studied analytical vectorial structure of Laguerre- 
Gaussian beam in the far field firstly [9] . Afterwards, Deng and Guo explored analytical 



vectorial structure of radially polarized light beams 



lOj. In 2008, Wu et al and Zhou et al 



investigated vectorial structure of hollow Gaussian beam almost in the same time ll|, Il2 |. 



In fact, much work has been done with the vectorial structure of all kinds of beams in the 
far field Q- Q. 

It is well known that the paraxial approximation is no longer valid for beams with a 
large divergent angle or, especially, a small beam spot size that is comparable with the light 
wavelength [l^]. Therefore, rigorous non-paraxial and vectorial treatments are necessary. We 
can approach non-paraxial propagation of beams in terms of vector angular spectrum method 
of electromagnetic field. According to vector angular spectrum method of electromagnetic 
field, the general solution of the Maxwell's equations is composed of the transverse electric 
(TE) term and the transverse magnetic (TM) term. In the far field, the TE and TM terms 
are orthogonal to each other and can be detached. 

To the best of our knowledge, the research on the vectorial structure of non-paraxial 
four-petal Gaussian beam in the far field based on vector angular spectrum method and sta- 
tionary phase method has not been reported elsewhere. In this paper, the far-field vectorial 
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properties of non-paraxial four-petal Gaussian beam have been studied by means of vector 
angular spectrum method and stationary phase method. Based on the analytical vectorial 
structure of the FPGBs, the energy flux distributions of TE term, TM term and the whole 
FPGBs are also investigated, respectively. 



II. ANALYTICAL VECTORIAL STRUCTURE IN THE FAR FIELD 



Let us consider a half space z > filled with a linear homogeneous, isotropic, noncon- 
ducting medium characterized by electric permittivity e and magnetic permeability fi. All 
the sources only lie in the domain z < 0. The electric field distribution is known at the 
boundary plane z = 0. For convenience of discussion, we consider a non-paraxial FPGBs 
with polarization in x direction, which propagates toward the half space z > along the z 



axis. The initia 
be written by [l| 



transverse electric field distribution of the FPGBs at the z = plane can 



\ 2n / 2 I 2 

xy \ ( X + y 



Ex{x,y,0) = Gn { — ] exp —] ,n = 1,2,3 ... , (1) 

Ey{x,y,0) = 0, (2) 

where n is the beam order of the FPGBs; G„ is a normalized amplitude constant associated 
with the order of n; wq is the 1/e^ intensity waist radius of the Gaussian term. The time 
factor exp{—iut) has been omitted in the field expression. Fig.[T]shows intensity distributions 
of four-petal Gaussian beams at the initial plane z = for n = 1,5,9 and 13, respectively. 
From Fig. [H it can be seen that the intensity distributions is composed of four equal petals. 
The distance of small petals increases when beam order increases. Here, Wq is taken by 
A in the calculation, which is comparable with the wavelength. So, this is a non-paraxial 
problem. 

In terms of Fourier transform, the vectorial angular spectrum of electric field at the z = 



plane is expressed as [ll| 

Ax{p,q) = -^j j E^{x,y,0)exp[-ik{px + qy)]dxdy, (3) 

^yiv^Q) = -j2 j j Ey{x,y,0)eyip[-ik{px + qy)]dxdy, (4) 
where A denotes the wave length in the medium related wave number by A; = 27r/A. Substi- 
tuting Eqs. and (E]) into Eqs. (E]) and (HD, we find that 





FIG. 1: (Color online) Normalized intensity distributions of FPGBs with different beam order n at 
z = plane based Eq. ([1]). (a) n = 1, (b) n = 5, (c) n = 9, (d) n = 13. 
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(5) 



where iFi{-; ■) denotes confluent hypergeometric function, r(-) denotes the Gamma func- 
tion. It is well known that Maxwell's equations can be separated into transverse and longi- 
tudinal field equations and an arbitrary polarized electromagnetic beam, which is expressed 
in terms of vector angular spectrum, is composed of the transverse electric (TE) term and 
the transverse magnetic (TM)term, namely. 



(6) 



H{r) = HTE{^ + HTM{^ 



(7) 



where 



ETE{r) = I I -M3;(p,g)](ge^ -pcy) 

J J — oo V ' Q 

X exp{iku)dpdq, (8) 



oo 

\ 1^ J J -oo P ~r Q 

X exp{iku)dpdq, (9) 



and 
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y j-ooV +q 7 

X exp{iku)dpdq, (10) 



HTuir) 




[pA^{p, q) + gv43^(p, g)]7^(ge^ - pe^, 
) " 7 



X exp{iku)dpdq. 



(11) 



r = xCa; + + 2:62 is the displacement vector and e^, e^, denote unit vectors in the x, y, 
z directions, respectively; m = pa; + + 72; 6^ = + g^; 7 = a/1 — — g^, if + g^ < 1 
or 7 = iy^p^ + g^ — 1, if p^ + g^ > 1. The value of p^ + g^ > 1 corresponds to the evanescent 



wave which propagates along the boundary plane but decays exponentially along the positive 
z direction. 

In the far field framework, the condition k{x'^ +y'^ + z'^Y^'^ — )• 00 is satisfied due to z is big 
enough. Moreover the contribution of the evanescent wave to the far field can be ignored. 
By virtue of the method of stationary phase 

fly 

, the TE mode and the TM mode of 

the electromagnetic field can be given by 
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r = A/a;2 + + 2;^. Eqs. (fT2l) - (|T5l) are analytical vectorial expres- 
sions for the TE and TM terms in the far field and constitute the basic results in this paper. 
From Eqs. ([IS]), we find that 



where p = ^J x^ ^ y 



2 . 



ETi,(f)-E™(r) =0, (16) 

Ete{^ ■ HTM{r) = 0. (17) 

According to Eqs. (fT6l) and (fT7|l . the TE and TM terms of FPGBs are orthogonal to each 
other in the far field. 



III. ENERGY FLUX DISTRIBUTIONS IN THE FAR FIELD 

The energy flux distributions of the TE and TM terms at the z = const plane are 
expressed in terms of the z component of their time-average Poynting vector as 

{S,)te = \Re[E*TE X Hte].. (18) 

{S,)tm = \Re[E*TM X HtmU (19) 

where the Re denotes real part, and the asterisk denotes complex conjugation. The whole 
energy flux distribution of the beam is a sum of the energy flux of TE mode and TM mode, 
namely, 

(5,) = + (5,)tm, (20) 
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Substituting Eqs. 1^- (IIS]) into Eqs. dH])- (HH) yields 
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Therefore, the whole energy flux distribution of FPGBs in the far field is given by 
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FIG. 2: (Color online)Normalized energy flux distributions of FPGBs at the plane z = 500A for 
beam order n = 1. (a) The TE term, (b) The TM term, (c) The whole beam. 

The normalized energy flux distributions of the TE term, the TM term and the FPGBs 
at the plane z = 500A for different beam order n = 1,5,9,13 versus x/X and y/X are 
illustrated by Figs. [2]- El The used parameter is t^o = A. As can be seen from Figs. [2]- El 
the FPGBs split into a number of small petals in the far field, which differs from its initial 
four-petal shape. The FPGBs with beam order n is not a pure mode, which can be regarded 
as a superposition of two dimensional Hermite- Gaussian modes and different modes 
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FIG. 3: (Color online) Normalized energy flux distributions of FPGBs at the plane z = 500A for 
beam order n = 5. (a) The TE term, (b) The TM term, (c) The whole beam. 





FIG. 4: (Color online)Normalized energy flux distributions of FPGBs at the plane z = 500A for 
beam order n = 9. (a) The TE term, (b) The TM term, (c) The whole beam. 

evolve differently within the same propagation distance. The overlap and interference in 
propagation between different modes result in the propagation properties of the FPGBs in 
the far field. Furthermore, the number of petals is determined by the parameter n. The 
number of petals in the far field gradually increases when the parameter n increases, which 
has potential applications in micro-optics and beam splitting techniques, etc. Note that 
diameter of the central beam spot decreases when beam order n increases. This phenomenon 
has been discussed in some previous researches 



IV. CONCLUSIONS 



In summary, the vectorial structure of the non-paraxial four-petal Gaussian beam in the 
far field is expressed in the analytical form by using the vector angular spectrum method 
and the stationary phase method. The electric field and the magnetic field of the four-petal 
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FIG. 5: (Color online)Normalized energy flux distributions of FPGBs at the plane z = 500A for 
beam order n = 13. (a) The TE term,(b) The TM term, (c) The whole beam. 

Gaussian beam is decomposed into two mutually orthogonal terms, i.e., TE term and TM 
term. Based on the analytical vectorial structure of FPGBs, the energy flux distributions of 
the TE term, the TM term and the whole beam of FPGBs are derived in the far-field and 
are illustrated by numerical examples. The number of petals and diameter of the central 
beam spot in the far field are determined by the beam order n. The potential applications of 
the FPGBs are deserved investigation. This work is important to understand the theoretical 
aspects of vector FPGBs propagation. 
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